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ABSTRACT 

The assumption of self-preservation allows for an analytical determination of the 
energy decay in isotropic turbulence. Batchelor (1948), who was the first to carry out a 
detailed study of this problem, based his analysis on the assumption that the Loitsianskii 
integral is a dynamic invariant — a widely accepted hypothesis that was later discovered 
to be invalid. Nonetheless, it appears that the self-preserving isotropic decay problem 
has never been reinvestigated in depth subsequent to this earlier work. In the present 
paper such an analysis is carried out, yielding a much more complete picture of self- 
preserving isotropic turbulence. It is proven rigorously that complete self-preserving 
isotropic turbulence admits two general types of asymptotic solutions: one where the 
turbulent kinetic energy K ~ t ~ 1 and one where K ~ t a with an exponent ot > 1 that is 
determined explicitly by the initial conditions. By a fixed point analysis and numerical 
integration of the exact one-point equations, it is demonstrated that the K ~ t power 
law decay is the asymptotically consistent high-Reynolds-number solution; the K ~ t 
decay law is only achieved in the limit as t — ► oo and the turbulence Reynolds number 
vanishes. Arguments are provided which indicate that a/(~ f -1 power law decay is the 
asymptotic state toward which a complete self-preserving isotropic turbulence is driven at 
high Reynolds numbers in order to resolve the imbalance between vortex stretching and 
viscous diffusion. Unlike in previous studies, the asymptotic approach to a self-preserving 
state is investigated which uncovers some surprising results. 


•Research supported by the National Aeronautics and Space Administration under NASA Contract 
No. NAS1-18605 while the author was in residence at the Institute for Computer Applications in Science 
and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665-5225. 
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1. INTRODUCTION 


Despite the fact that isotropic turbulence constitutes the simplest type of turbulent 
flow, it is still not possible to render the problem analytically tractable without the 
introduction of additional hypotheses. The idealization of self-preservation — wherein 
the two-point double and triple longitudinal velocity correlations are assumed to admit 
self-similar solutions with respect to a single length scale L(t) — has served as a useful 
hypothesis since its introduction by von Karman and Howarth (1938). In another classic 
paper that followed, Batchelor (1948) studied the energy decay in self- preserving isotropic 
turbulence in considerable detail. He concluded that the only complete self-preserving so- 
lutions that were internally consistent existed at low turbulence Reynolds numbers where 
the turbulent kinetic energy K ~ t ~ 5/ 2 - a power law consistent with the final period 
of decay. Batchelor (1948) also found a self-preserving solution to the Karman-Howarth 
equation in the limit of infinite Reynolds numbers for which Loitsianskii’s integral was an 
invariant. This solution - wherein L(t) is the integral length scale A and K ~ t~ w/1 - 
was put forth by Batchelor as the only complete self-preserving solution at high Reynolds 
numbers. Of course, additional partial self- preserving solutions were shown by Batchelor 

to exist in other Reynolds number regimes. 

Objections were later raised against the use of the Loitsianskii integral as a dynamic 
invariant: at high Reynolds numbers this integral can be shown to be a weak function of 
time (see Proudman and Reid 1954 and Batchelor and Proudman 1956). Saffman (1967) 
proposed an alternative dynamic invariant which yields a K ~ U 6/5 power law decay in 
the infinite Reynolds number limit (see Hinze 1975). While the results of Batchelor and 
Saffman formally constitute complete self-preserving solutions to the mviscid Karman- 
Howarth equation, it must be kept in mind that they only exhibit partial self-preservation 
with respect to the full viscous equation. Namely, there is self-preservation only for the 
range of energy containing eddies with integral length scale A (here, A/ A << 1 where A is 
the Taylor microscale). These two solutions have been widely accepted in the turbulence 
literature as the predicted decay laws for self-preserving isotropic turbulence at high 
Reynolds numbers. 

Implicit in the analysis of Batchelor (1948) is the existence of a complete self- preserving 
solution consistent with high Reynolds numbers - namely a I< ~ t~ l power law decay. 
The collapsing length scale L(t) for this full self-preserving solution is necessarily the 
Taylor microscale (i.e., for any complete self-preserving solution of the viscous Karman- 
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Howarth equation we must have L oc A). This solution — which was postulated earlier by 
Dryden (1943) — was dismissed by Batchelor on the grounds that Loitsianskii’s integral 
was not a dynamic invariant. While this result by Dryden has been mentioned subse- 
quently in the literature (c.f. Hinze 1975, Monin and Yaglom 1975 and Korneyev and 
Sedov 1976), it has largely been disregarded by the turbulence community. The reason 
for this appears to be two-fold: (a) a K ~ i -1 power law decay has not been observed 
in the most accurate isotropic decay experiments, and (b) since A/ A — ► 0 as Re — > oo, 
doubts can be cast on the suitability of the Taylor microscale as the collapsing length 
scale of the energy containing eddies. 

Recently, George (1987, 1989) revived this issue concerning the existence of complete 
self-preserving solutions in isotropic turbulence. In an interesting paper he claimed to 
find a complete self-preserving solution, valid for all Reynolds numbers, in which the 
kinetic energy decayed as K ~ t~ a with a determined by the initial conditions. George 
(1987) — who based his analysis on the dynamic equation for the energy spectrum rather 
than on the Karman-Howarth equation — made no explicit mention of the complete self- 
preserving I\ ~ solution. Strictly speaking, the solution presented by George was an 
alternative self-preserving solution to that of Karman and Howarth (1938) and Batchelor 
(1948) since he relaxed the constraint that the triple longitudinal velocity correlation be 
self-similar in the classical sense. 

The purpose of the present paper is to address the issue of complete self-preservation 
in an effort to clarify the following basic questions: 

(1) What is the complete self-preserving solution for isotropic turbulence at high 
Reynolds numbers, 

(2) What detailed predictions does this solution yield for the energy decay, particularly 
during the initial approach to a state of complete self-preservation, and 

(3) Is this solution compatible with the results of physical experiments and alternative 
theoretical approaches? 

In so far as the first two points are concerned, it will be shown unequivocally that the 
only complete self-preserving solution that is consistent with a high-Reynolds-number 
isotropic turbulence has a K ~ t -1 asymptotic power law decay. Unlike previous studies, 
this is demonstrated in a straightforward manner based on a fixed point analysis of the 
one-point equations. This analysis leads to an interesting interpretation of the physical 
significance of a K ~ < _1 power law decay and allows us to examine small departures 
from a state of complete self-preservation. The detailed predictions of this complete self- 
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preserving solution - which, to the best of our knowledge have never been examined in 
depth in the literature — will be compared with the results of physical experiments and 
alternative theoretical approaches in the sections to follow. 


2. THEORETICAL BACKGROUND 


We will consider isotropic turbulence governed by the incompressible Navier-Stokes 
equations 

"" "" ( 1 ) 


dui dui dp 2 

— 1 + Uj-r — = — 5 + I/V Ui 
dt 3 dxj dxi 


OXi 


( 2 ) 


where Uj is the velocity vector, p is the pressure and v is the kinematic viscosity. The two 
point double and triple longitudinal velocity correlations, denoted by f(r,t) and k(r,t), 
respectively, are defined in the standard way: 


„ .X ti(x,t)u(x+r,i) 

/(’-.o = p 


( 3 ) 


*( r »0 = 


u 2 (x,t)u(x + r, t) 


- , ( 4 ) 

(it 2 ) 3 /2 

where u is any component of the velocity, x and x + r are any two spatial points separated 
by a distance r =| r | in the direction of u, and an overbar denotes a spatial average (c.f. 
Hinze 1975). For isotropic turbulence, / and k satisfy the Karman-Howarth equation 

d&f) _ f~2 


dt 




( 5 ) 


which is obtained directly from the Navier-Stokes equations. The turbulent kinetic energy 
ft — 1 uftTi is a solution of the differential equation 


K = -< 


161 


where 

t = v = v uj 2 


( 7 ) 


is the turbulent dissipaton rate, u>< is the vorticity vector and a ? is the enstrophy. The 
turbulent dissipation rate is a solution of the differential equation 
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where 


Sk 



-A 3 


d 3 k 


dr 3 


r=0 


(9) 


G=X 4 



( 10 ) 


D K 2 . / \QuK\ lt 2 

Rt ~ vt ’ A ~ ( c ) ( n ) 

are the velocity derivative skewness, the coefficient of the enstrophy destruction, the 
turbulence Reynolds number, and the Taylor microscale, respectively. Equations (6) and 
(8) — which are rearranged versions of those that appear in Karman and Howarth (1938) 
and Batchelor (1948) — are obtained, respectively, by setting r = 0 in (5) and in the 
equation obtained by differentiating (5) twice with respect to r. Since Sk = SW(<) and 
G = G(t) are directly related to the correlations / and k (which cannot both be obtained 
from the Karman-Howarth equation (5)) it is clear that the problem of isotropic decay 
is not closed. In order to achieve closure, additional assumptions must be made such as 
the one of complete self-preservation that will be discussed in this paper. 

For an isotropic turbulence to be self-preserving in the sense of Karman and Howarth 
(1938) and Batchelor (1948), we must have 


/(’■>*) = /(r/I) 

(12) 

ii 

(13) 


where L = L(t) is a uniquely specified similarity length scale. For it to exhibit complete 
self-preservation, all scales of the turbulence — namely, the full range of 0 < r < oo 
— must decay according to (12) and (13) (partial self-preservation is satisfied if (12) - 
(13) only apply to some restricted range of 0 < r < r max ). We will focus our attention 
on complete self-preserving solutions in the analysis to follow. In view of the identity 
(Batchelor 1948) 


it follows from (11) that 


e = -10 vK 

a 2 m 

dr 2 


d 2 f 
dr 2 n 

J r=0 

(14) 

= -l. 

(15) 


J r=0 
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( 16 ) 


Hence, for any complete self-preserving isotropic turbulence, we must have 

^/"( 0 ) = -1 

from which it can be concluded that 


L oc A (17) 

since /"( 0) is a constant. It therefore follows that the Taylor microscale is the only 
similarity length scale that can yield complete self-preserving solutions to the full viscous 
equations of motion for isotropic turbulence. 

Without any loss of generality for a complete self-preserving isotropic turbulence we 
may set L = A and then substitute (12) - (13) into (9) - (10), respectively, to get 

Sk = —k'"( 0) = constant 

G = f tv ( 0) = constant 

where a prime denotes a derivative with respect to rj = r/A. Consequently, 


Sk = S Ko , G = G 0 


(18) 


where the notation (•)<, denotes the initial value. The substitution of (18) into (6) and 
(8) then yields the transport equations 


K = 


e = 


m s ^ e K~h a ^ 


(19) 

( 20 ) 


for complete self-preserving isotropic turbulence. This is a closed system for the deter- 
mination of I{ and e once K 0 ,c 0 ,Sk o and G 0 are provided. To simplify the subsequent 
presentation, the quantity G 0 - which is the coefficient of the destruction of enstrophy 
term in (20) - will henceforth be referred to as the “palinstrophy coefficient” following 
the terminology used by Lesieur (1990). Accordingly, the assumption of complete self- 
preservation is seen to lead to closure in the following sense: if initial conditions for the 
skewness and the palinstrophy coefficient are provided — in addition to initial conditions 
for K and e — then the energy decay can be calculated explicitly for all later times. 

For complete self-preserving isotropic turbulence, the Karman-Howarth equation (5) 
takes the form 


Kf 


X df h K _ 4 d(rj 4 k) K _ 4 d ( 4 dp 


0 


( 21 ) 
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or, equivalently, 


10/ + I T] 




dr] \ dr] 


, tin 

+ %{3 G ° 


i) = Rx (gSjf, 


df_ 4 djrfjtf 

V dr] V dr] 


( 22 ) 


after replacing A using (11), (19) and (20) with R\ = {u 2 ) l l 2 \/u = ^20/3 R^ 2 . Equation 
(22) will have a solution if R\ — constant as first noticed by Dry den (1943); this is a 
K ~ t~ x power law decay. However, (22) also has solutions where R\ = R\{t) when 
separability is invoked. The separability condition implies that each side of (22) is equal 
to zero individually, yielding differential equations from which explicit solutions for / 
and k may be determined depending on the choice of Sk 0 and G 0 . These solutions were 
first discovered by Sedov (1944) and later compared with experimental data by Korneyev 
and Sedov (1976). The particular case for which G a = 3 so that / is Gaussian — which 
formally corresponds to the final period of decay — was considered in detail by Bernard 
(1985). We will briefly examine (22) later to establish the consistency of the present 
results with those of previous studies. However, our analyses will be based on a fixed 
point analysis and direct numerical integration of (19) - (20). This will allow us to 
consider small departures from a self-preserving state as will be demonstrated later. 


3. FIXED POINT ANALYSIS AND NUMERICAL RESULTS 


In order to carry out a fixed point analysis of (19) and (20), we will combine these 
equations into a single transport equation for the turbulence Reynolds number R t . Since, 


it follows that 


Rt = - 


• 2 KI< K 2 . 

Rt = 

ve vt l 


2K 7 c rw K 4- 7 r- K 

--jjE SK °'J Rt V + Ti G ‘V 


(23) 


(24) 


If the transformed dimensionless time r — defined by the relation dr = ( e/K)dt — is 
introduced into (24), we obtain the equation 


f = R -{Ti G °- 2 -m S ^ R ‘ ] - 


(25) 


The fixed points of (25) are obtained by setting dR t /dr = 0 which yields the equation: 

R ‘"{h G °- 2 -m SK -J* :: ) =0 (26) 
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where (•)<» denotes the equilibrium value in the limit as r — > oo. Equation (26) has the 
solutions: 

Rt x = 0 (27) 

for yeG 0 < 2, and 


ft. = -“t 

\ 57Ts 


■hGo — 2 


Sk 0 


(28) 


for y^G 0 > 2. It is a simple matter to show that the fixed points (27) - (28) are stable 


nodes that attract all initial conditions K a and e Q . It is also evident from (28) that in order 
to have an equilibrium high-Reynolds-number isotropic flow field^ under self-preserving 
conditions it is necessary that G 0 ~ \fRtZ- 

The fixed point R too — 0 is associated with asymptotic solutions of K and e that 
satisfy the differential equations 

K = -e (29) 


7 e 2 

* = -ii a °T 


Equations (29) - (30) have the exact solution 


*-*-"*S) 


-K0 


—o — l 


(30) 

(31) 

(32) 


where a = 1/(ttG 0 — 1) > 1. For large eddy turnover times t 0 tj K 0 1, we have 


I< ~ t~ a (33) 

Rt ~ r Q+1 (34) 

(and, hence, K goes to zero faster than Rt), so it follows that (31) - (32) can only be 
asymptotic solutions in the limit as K — *■ 0 and t — *• oo. Therefore, within the framework 
of complete self-preservation, an asymptotic power law decay where K ~ t~ a (with 
a > 1) is associated with the fixed point R too = 0 and is only formally consistent with 
the final period of decay. This is borne out in the subsequent computations. 

During the final period of decay there is considerable evidence (Batchelor and Townsend 
1948a) indicating that /(t/) = exp (— r} 2 /2) (i.e., a Gaussian) in which case (10) implies 

tBy a high Reynolds number isotropic turbulence we mean the case where R t '$> 1; for a low Reynolds 
number isotropic turbulence, Rt — 0(1). 
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that G 0 — 3 and, consequently, that a = 5/2. The same result is also reached by 
assuming constancy of the Loitsianskii integral 

yoo 

u 2 J r 4 f(r,t)dr = constant (35) 

during decay, which appears to be an acceptable assumption for the final period. In 
particular, from (12) it follows that (35) is equivalent to (c.f. Hinze 1975) 

roo „ 

u 2 X 5 J t} 4 f (rj)drj = constant (36) 

so that (11) and (36) imply that 

R7/2 

—jjY — constant. (37) 

When (37) is combined with (31) and (32) it follows that a — 5/2 — the celebrated 
Batchelor (1948) result. 

Now we will show that the non-zero fixed point (28) is consistent with high-Reynolds- 
number isotropic turbulence. The substitution of (28) into (19) - (20) yields the equations 


<o 

1 

II 

(38) 

-4 

which have the exact solution 

(39) 


(40) 

--('*gr 

(41) 

(namely, a t 1 power law decay for the turbulent kinetic energy), 
that (see Batchelor and Townsend 1948b) 

Since it can be shown 

n 30l^ 1 
° ~ 7 e/K 

(42) 

it follows that G is a ratio of turbulent to dissipative time scales, 
to show that (see Hinze 1975) 

It is a simple matter 

_ Vf 0 °° K 4 E(K,t)dK 
UkUk fo°° K 2 E(n,t)dK 

(43) 
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where E(n,t) is the energy spectrum and k is the wavenumber. For complete self- 
preservation, E(K,t) scales with the Taylor microscale, from which it follows that 


and, hence, that 


du/j duSj 

” dxj dxj V £ 

— rsj 

UJkUJJi A 2 I< 


G ~ constant 


(44) 

(45) 


as shown earlier. However, (43) is a correlation dominated by the high wavenumbers 
(i.e., small scales) and it would therefore seem more reasonable that E(k, t ) should scale 
with the Kolmogorov length scale Ik = If this is the case, then 


dujj duti 

dxj dxj V 

- ■■■ ■ 

^k^k 

and 

G~~~y/Rf (47) 

l K € 

Equations (45) and (47) appear to be contradictory; however, they are not, since in a 
K ~ t~ 1 power law decay 

Rt — constant. (48) 

Furthermore, R t = constant resolves the imbalance between the two terms on the right- 
hand side of the dissipation rate equation (20) since the first term (i.e., the vortex stretch- 
ing) is initially of 0(\fRt) while the second term (i.e., the viscous destruction) is of 0(1). 
This leads us to the following physical interpretation: A K ~ t~ l power law decay is 
the asymptotic state toward which a self-preserving isotropic turbulence is driven at high 
Reynolds numbers in order to resolve the fundamental imbalance between vortex stretching 
and viscous diffusion. In the process of resolving this imbalance, compatibility with Kol- 
mogorov scaling is achieved for the small scale correlations. Since this consistency with 
Kolmogorov scaling — which, on physical grounds should be satisfied at high Reynolds 
numbers — is achieved when ~ > 2, it is clear that the non-zero fixed point 

(28) is the asymptotic solution for high-Reynolds-number self-preserving flows. 

We will now examine numerical solutions of (19) - (20) for complete self-preservation. 
With the possible exception of the recent work of Bernard (1985), we have not seen 
detailed numerical results published on the decay of K and e in self-preserving isotropic 
turbulence. An examination of these results will amplify the points discussed in this 
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section and will shed some interesting new light on how the self-preservation assumption 
compares with experiments. In Figures 1(a) - (c), the decay of the turbulent kinetic 
energy is shown for three initial turbulence Reynolds numbers ( R to = 1000, 5000 and 
10,000) and three different initial conditions for G (i.e., C t2 = 1.92,5.0, and 8.0 where 
= yeG 0 ). For these calculations, as well as the ones to follow, K* = K/K a ,t* = 
c 0 t/K 0 and the skewness 

Sk 0 — 0.5 (49) 

which is in close proximity to the values obtained from physical experiments in this range 
of Reynolds numbers. From these figures it is clear that the self- preserving solution has 
an initial transient where the kinetic energy is fairly flat; then the kinetic energy begins 
to asymptote from above to a power law decay as evidenced by a straight line on these 
logarithmic plots. Two conclusions can be drawn from these results. First, for C t2 < 2 
and R to > 1, the kinetic energy does not asymptote to a t~ a power law decay until 
after an extremely large number of eddy turnover times by which time the turbulence 
has decayed to a tiny fraction of its initial intensity. Second, for C t2 > 2, the kinetic 
energy asymptotes to a t ' power law decay within a few eddy turnover times; however 
unless C t2 ~ \fTZT 0 (consistent with Kolmogorov scaling), the turbulence intensity will 
drop precipitously before this asymptotic state is achieved. 

To further illustrate these points, the computed turbulent kinetic energy is compared 
with its corresponding asymptotic power law decay for increasing values of C C2 in Figures 
2(a) - (c) which are for initial turbulence Reynolds numbers of 1000, 10,000 and 100,000, 
respectively. It is clear from these figures that for C C2 < 2, the turbulent kinetic energy 
does not reach its asymptotic power law decay even after 100 eddy turnover times! How- 
ever, for C t2 sufficiently larger than 2, the turbulent kinetic energy asymptotes to a t~ l 
power law decay within a few eddy turnover times. This can be seen even more vividly in 
Figure 3 where a* = d(log K*)/d{\og t*) is plotted as a function of log t* for R to = 1000. 
If there is an asymptotic power law decay, this derivative will asymptote to the exponent 
of the decay law. It is clear that for C C2 = 7, an exponent of 1 is approached quickly; 
however, for C £2 = 1.83 and 1.92 (initial conditions which ultimately yield a power law 
decay with an exponent of approximately 1.2 and 1.1, respectively) an asymptotic state 
is not achieved even after 100 eddy turnover times. Furthermore, for C t2 < 2 and large 
initial turbulence Reynolds numbers R to 1, there is a precipitous drop in the turbulent 
kinetic energy before a power law decay is achieved; this is due to the early transient 
when vortex stretching causes a considerable rise in the dissipation (see Figures 4 (a) - 
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(b))- 

Since the self-preserving solutions for C (2 < 2 only asymptote to a power law decay 
in the limit as K — * 0 and t —* oo, it is reasonable to associate them exclusively with 
the final period of decay. Experiments tend to indicate that the final period of decay is 
entered for R t < 1 wherein the exponent of the decay is approximately 2.5 (c.f. Hinze 
1975). As noted earlier, this decay law is obtained asymptotically for self-preserving 
isotropic turbulence if C (2 = 1.4 — a result obtained by invoking Loitsianskii s invariant. 
In Figure 5, the decay of the turbulent kinetic energy when C t2 = 1.4 is shown for the 
initial turbulence Reynolds numbers Rt a — 0.25 and 1.0. It is clear from these results 
that for R to < 1 the solutions begin to asymptote to the power law decay (31) with an 
exponent a = 2.5. Since C t2 must equal 1.4 during the final period of decay — and, since 
C t2 is a constant for self-preserving isotropic turbulence — it follows that the entire decay 
process from high Reynolds number initial conditions to the final period of decay cannot 
be described within the framework of complete self-preservation. This conclusion results 
from the fact that the only consistent self-preserving solution at high initial Reynolds 
numbers yields a K ~ t~ * power law decay wherein Rt asymptotes to a constant a 
state of affairs that precludes the description of the later stages of decay. Furthermore, 
the value of C t2 = 1.4, which describes the final period of decay, yields unphysical results 
for the early stages of a high Reynolds number isotropic turbulence (i.e., it predicts an 
early time transient where there is a precipitous drop in the turbulent kinetic energy; 
see Figure 6). In order to describe the entire decay process of a high Reynolds number 
isotropic turbulence, G as well as Sk must vary with time a possibility that is precluded 
by the assumption of self-preservation which renders them to be constant. 

4. COMPARISONS WITH ALTERNATIVE THEORETICAL ANALYSES 
AND EXPERIMENTS 

The results derived in the previous section are consistent with those of Batchelor 
(1948) for low turbulence Reynolds numbers; however, our high Reynolds number asymp- 
totic solution yields I\ ~ t -1 whereas in Batchelor’s solution A ~ t 10 / 7 . The reason 
for this difference is simple: as alluded to earlier, Batchelor also found the K ~ t _1 
solution but dismissed it as a viable result since Loitsianskii’s integral was not a dy- 
namic invariant. Interestingly enough, an earlier experimental study by Batchelor and 
Townsend (1948b) yielded results that were far more suggestive of a K ~ t~ l rather 
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than a I< ~ r 10/7 power law decay. Despite the fact that Batchelor (1948) states that 
the K ~ t ~ 10 / 7 decay law is a complete self-preserving solution, in reality it is only a 
partial self-preserving solution corresponding to the inviscid Karman-Howarth equation 
(see Hinze 1975). It is our view that since the K ~ t~ l asymptotic decay law is a formal 
solution to the full Karman-Howarth equation, it should not be dismissed unless it is in 
incontrovertible contradiction of experiments or other exact theoretical results. 

As mentioned earlier, Dryden (1943) postulated a I< ~ t~ l power law decay based 
on a direct analysis of the Karman-Howarth equation. He observed — as is evident from 
(22) — that the Karman-Howarth equation will allow for self-similar solutions if 


R\ = constant 


which yields K ~ t 1 as a direct consequence of (6). However, there are other temporally 
varying solutions to (22); complete self-preservation only requires that R x asymptote to 
a constant. Sedov (1944) studied solutions of (22) obtained by applying the separability 
constraint 


10 / + 2t]' 




(7 df _ 4 d(rj 4 k)\ 

7 drj j~° ( 52 ) 

which renders R\ = Rx(t) consistent with (19) - (20). Solutions to (51) - (52) have 
not been studied in great depth subsequent to Sedov (1944) who showed that /(?/) = 
w h ere AT is the confluent hypergeometric function and 7 = 7q q _ |_ 
Batchelor (1948) expressed concern over the fact that this solution leads to a unique 
determination of both / and k ; however, although he suspected that the Sedov solution 
was unphysical, he stated that he was “not able to find any definite anomalies”. In 
the limit as R\ — * 0 it is clear that (51) is a direct consequence of the Karman-Howarth 
equation. Consequently, it is not surprising that the Sedov solution for the final period of 
decay yields physically interesting solutions as recently demonstrated by Bernard (1985). 
However, it will now be shown definitively that the Sedov solution yields unphysical 
results at high turbulence Reynolds numbers. In figure 7, the results of a numerical 
solution of (51) for / are shown for a variety of values of G 0 ranging from 3 to 60. For 
G 0 3 it can be shown analytically that f = exp( — ?^ 2 /2) yielding an energy spectrum 
of the form 

E *( k *) = exp(-AC* 2 /2) (53) 
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( 54 ) 


where k* = k\ and E* = E/u 2 A. Equation (53) is obtained from the identity 

1 roo „ 

E*(k*) = — / f{v)( K *V s ^ n K *V ~ t/ 2 cos K*7])dr] 

(c.f. Batchelor 1953). This result — which has E*(k*) ~ ft* 4 at low wavenumbers and 
has 2£*(ft*) decaying exponentially at high wavenumbers — is consistent with established 
results on the final period of decay (c.f. Hinze 1975). However, for sufficiently large (7 0 , 
it is a simple matter to show from (51) that 


f(v) ~ V 7 


(55) 


for rj 1. This explains why f(t}) is so slow to asymptote to zero when G 0 > 10 in 
figure 7. In fact for 

—G 0 > 2 

it follows from (54) and (55) that E*(k*) becomes singular. From (28) it can then be 
concluded that the Sedov solution will yield a singular energy spectrum when 

Rtco > 0 

It is thus clear that the Sedov solution yields untenable results, at high Reynolds numbers, 
for the double and triple two-point longitudinal velocity correlations. 

The major deficiency with the approaches of Dryden and Batchelor — as well as that 
of Sedov — lies in the use of the self-similar Karman-Howarth equation (22) which does 
not allow for the treatment of small departures from a state of complete self-preservation. 
Such small departures can be characterized by the perturbations 


Sk = S Ko + SS K (t) 

G = G 0 + 6G(t) 

/G 0 < 1 and (55/c(t), SG(t) 


(56) 


(57) 

where || SSk || /Sk 0 <. 1, || $G || /G 0 <C 1 and (55/c(t), 6G(t) — > 0 as t — » oo. The 
substitution of (56) - (57) into (6) and (8) yields the governing equations for small 
departures from a state of self-preservation. If we denote SK and St as the departures 
from the self-preserving solutions K and t obtained from (19) - (20), it follows that for 
the perturbations (56) - (57) we will have 


due to the fact that (19) - (20) have fixed points that are stable nodes (c.f, Guckenheimer 
and Holmes 1986). Consequently, (19) - (20) will yield an excellent approximation for 
isotropic decay when there are extremely small departures from a state of complete self- 
preservation. In contrast to this nice behavior, the Karman-Howarth equation becomes 
indeterminate when subjected to infinitesimal perturbations from a self-preserving state. 
Hence, it appears that the one-point equations (19) - (20) form a broader basis for the 
analysis of the energy decay of self-preserving isotropic decay than does (22). 

The general solution to the complete self- preserving isotropic decay equations (19) - 
(20) at high Reynolds numbers is shown schematically in figure 8 (this is for the physically 
significant case where R to > Rt^ so that the turbulence Reynolds number decays). There 
is an early time transient (region AB) where the turbulent kinetic energy is flat; it is 
eventually followed by the asymptotic region CD where K ~ t~ x . These two regions are 
connected by the overlap region BC . The initial transient AB evolves on the Kolmogorov 
time scale yjvje during which time there is a precipitous drop in the turbulence Reynolds 
number (see figure 9). On the other hand, the overlap region BC evolves on the turbu- 
lence time scale K /e; in this region the turbulence Reynolds number R t becomes close to 
R too , approaching it asymptotically from above. As a direct consequence of the perturba- 
tion analysis discussed above, the overlap region BC can be set into strong approximate 
agreement with the asymptotic approach to a state of complete self-preservation. These 
results have a direct bearing on how the complete self-preserving solution compares with 
physical experiments as we will soon see. 

It is widely believed that a K ~ t~ x asymptotic decay law is in violation of experi- 
mental data for isotropic turbulence. This experimental data (see Uberoi 1963, Kistler 
and Vrebalovich 1966, Comte-Bellot and Corrsin 1966, 1971, and Warhaft and Lumley 
1978) has yielded power law decays with exponents varying from 1 to 1.4 with a mean of 
approximately 1.25. However, great caution must be taken in using this data to dismiss 
the possibility of a K ~ t~ x asymptotic power law decay at high Reynolds numbers since 
most of this data is for a limited number of eddy turnover times (typically for e 0 tj K a < 4). 
If the self-preserving solution is examined for this same limited number of eddy turnover 
times it follows that the resulting solution can be fitted to an excellent degree of approx- 
imation by a power law decay with exponents in the range of 1 to 1.4 depending on the 
initial conditions; the lower the Reynolds number, the longer the solution takes to reach 
an asymptotic state and the larger the exponent is during the early stages of decay (see 
Figures 10 (a) - (c)). Consequently, if one examined in isolation the self-preserving solu- 
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tions for the first few eddy turnover times (with the short early time transient omitted), 
one could erroneously conclude that there was an asymptotic power law decay with an 
exponent in the range of 1 to 1.4 depending on the initial conditions; in reality, all of 
these solutions are asymptoting to a f -1 power law decay. The solutions shown in Figures 
10 (a) - (c) correspond to the overlap region BC shown in Figure 8 and, hence, can be 
associated with the asymptotic approach to a state of complete self-preservation. An 
argument has been raised recently by Walker and Corrsin (1985) and Walker (1986) that 
the physical experiments may not go far enough to see a t~ l power law decay. Unless 
the initial turbulence Reynolds number is extremely large, an asymptotic state may not 
be achieved in the first few eddy turnover times. In this regard it is interesting to note 
that the only extremely high- Reynolds- number experiment (i.e., Kistler and Vrebalovich 
1966) and large-eddy turnover time experiment (Walker 1986) did measure a K ~ t 1 
asymptotic power law decay. Consequently, existing experiments cannot rule out the 
possibility of a K ~ t~ l asymptotic power law decay at high Reynolds numbers and 
do not warrant the dismissal of the complete self- preserving solution discussed herein. 
Furthermore, Rosen (1981) recently derived a t -1 asymptotic power law decay based on 
alternative ideas from statistical mechanics. 

Now, we will address the interesting controversy generated recently by George (1987, 
1989). He claimed to find complete self-preserving solutions — with the Taylor microscale 
as the similarity length scale - that exist for all turbulence Reynolds numbers. These 
solutions were characterized by an asymptotic power law decay where the exponent is 
determined by the initial conditions. George arrived at this alternative self-preserving 
solution by relaxing the classical similarity constraint (13). He argued that the normal- 
ization of the two-point triple velocity correlation 

T(r,t) — u 2 (x, f)u(x + r,f) (58) 

by (u 2 ) 3 / 2 to form k(r, t) is arbitrary since its one-point contraction T(0,f) = u 3 is zero. 
(This stands in contrast to the formulation of /(r, t) which is obtained by normalizing the 
two-point double velocity correlation with its one-point contraction u 2 ). Consequently, 
George argued that constraint (13) should be replaced with the alternative constraint 

W = (59) 

where w(t) is a suitable weighting function. Then — from the definition of Sk in equation 
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(9) and the Karman-Howarth equation (22) — instead of the constraints 

Rx = constant , Sk = constant (60) 

which renders a f -1 power law decay, we get the constraints 

R\Sk = constant, R\wK ~ 3 ^ 2 — constant (61) 

which allows for the possibility of an alternative decay law within the general framework 
of complete self-preservation. From (61) it follows that 

S K °c R7 1/2 (62) 


and that 

w oc K 3/2 R; 1/2 , (63) 

which shows, incidentally, that w cannot be chosen arbitrarily. Since the proportionality 
constant in (62) must be Sk q R V* ~ and since G still remains a constant G 0 during the 
decay — this alternative self-preservation leads to the decay equations 


K = — e 


(64) 


( c e»!/2 € * n t 

1 = iTf!' *’*• K - I5 G “¥ 


(65) 


instead of (19) - (20). Equations (64) - (65) yield the closed form solution for the energy 
decay 

( 66 ) 


K = K 0 (l + ~ ot 


PK 0> 


where 

(67) 

This is indeed a power law decay with an exponent that depends on the initial conditions 
as claimed by George (1987, 1989). 

Although we do feel that George raised some interesting issues — and made a seem- 
ingly valid point concerning the arbitrariness of the normalization of the two-point triple 
velocity correlation — it is our opinion that there is a problem with this alternative 
self-preserving solution. By virtue of (61), George’s alternative scaling renders a zero 
skewness in the limit of infinite turbulence Reynolds numbers. Experiments indicate 
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that the skewness is an absolute constant of order one at extremely high Reynolds num- 
bers (see Van Atta and Antonia 1980). Furthermore, a direct physical consequence of a 
vanishing skewness Sk is the occurrence of zero transfer and, hence, no energy cascade. 
This can be seen easily from equation (59) which can be re-written as 

T{r,t)<xK 3 / 2 R; 1/2 T (£) 


which implies that the transfer T{r,t) — » 0 as Rt —+ oo. Vanishing transfer in the limit 
of infinite Reynolds numbers is an untenable physical result and hence this alternative 
self-preserving solution is not acceptable for high-Reynolds-number isotropic turbulence. 

Finally, a few comments are in order concerning the implications of these results 
for turbulence modeling. In the commonly used turbulence models, the dissipation rate 
equation is modeled as 

(«) 

for isotropic decay, where C e , is a constant (c.f. Launder and Spalding 1974, Speziale 
1991). Equation (68) is derived by invoking Kolmogorov scaling for G which requires 
that 

G = C\\J~Rt + C 2 (69) 


where C\ and C 2 &re constants; an equilibrium hypothesis is then made by which C\ 

so that the leading order part of the destruction of dissipation term annihilates 
the vortex stretching term in (8) yielding (68). In contrast to (68), the complete self- 
preserving solution has a dissipation rate transport equation of the general mathematical 
form 

‘e 4 
*K 


< Rt — — C t2 — 


e = C 


£3 


(70) 


where C C2 and C t3 are constants. Equation (70) can be also derived based on Kolmogorov 
scaling (69) when departures from equilibrium are allowed wherein C\ ^ ^ Sk • The 
addition of the unbalanced vortex stretching term in (70) allows for a better treatment 
of departures from equilibrium in several ways. First, as shown earlier in Figures 10 (a) - 
(c), the self-preservation model allows for the description of the initial stages of isotropic 
decay where the exponent of the decay law can vary mildly with the initial conditions 
a feature observed in physical experiments. On the other hand, the more commonly used 
model (68) predicts a universal decay law where K ~ t Ce 2 -1 for all Rt a and all t > 0. 
Second, the self- preserving solution can accommodate the limit of zero viscosity. In this 
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limit, it is a simple matter to show that (20) predicts a finite time enstrophy blow-up at 
the critical time 

6v/l5 
C 7 S Ko u>o 

where u 2 is the initial enstrophy. In Figure 11(a), the time evolution of the enstrophy 
corresponding to the complete self-preserving solution is shown for a variety of increasing 
R to ; it is clear that an enstrophy crisis is predicted for R to >■ 1 which eventually leads 
to a finite time enstrophy blow-up in the limit as v — » 0. These results are in excellent 
qualitative agreement with results obtained from EDQNM as illustrated in Figure 11(b) 
taken from Lesieur (1990). While the issue of a finite time enstrophy blow-up is still 
being debated by the turbulence community (c.f. Pumir and Siggia 1990), one thing is 
clear: the enstrophy grows dramatically when v — 0. In contrast to the results shown in 
Figures 11 (a) - (b), the commonly used dissipation rate model (68) erroneously predicts 
that the enstrophy is conserved in the inviscid limit, i.e. that 

u 2 = constant (72) 

when v = 0. It thus appears that the complete self-preserving solution allows for a 
better treatment of non-equilibrium isotropic turbulence that could be of future use in 
the development of improved turbulence models. 

5. CONCLUSIONS 

The energy decay for complete self-preserving isotropic turbulence has been re-examined 
from a basic theoretical and computational standpoint. Several interesting conclusions 
can be drawn from these results: 

(1) The nonlinear differential equations for the energy decay have two fixed points — 

R too = 0 and R toa = 0 — 2 ) 2 /Sj Co . The former fixed point is only achieved in the 

limit as t — > oo and hence is associated with the final period of decay. Consistent with 
the Batchelor (1948) result, a K ~ t~ 5 / 2 power law decay is obtained when Loitsianskii’s 
invariant or the Gaussianity of /(t/) is invoked. 

(2) The non-zero fixed point R toa — — 2) 2 /S] io is approached within a few 

eddy turnover times and gives rise to a K ~ f -1 asymptotic power law decay. It is the 
high-Reynolds-number asymptotic solution for a complete self-preserving isotropic tur- 
bulence. This solution appears to have been prematurely dismissed by Batchelor (1948) 
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purely on the grounds that Loitsianskii’s integral was not an invariant - a constraint 
which was later learned to be violated in isotropic turbulence when Rt a 1. 

(3) The structure of the high Reynolds-number self-preserving solution during the 
first few eddy turnover times was examined in detail. By a perturbation analysis, it was 
argued that these solutions can serve as an approximation for the asymptotic approach 
to a state of complete self-preservation. It was found that, depending on the initial 
conditions, the early time solutions could be fitted with a power law decay which has an 
exponent varying from 1.0 to 1.4 — a range of values that is compatible with existing 
experimental data. Consequently, existing experiments cannot rule out the possibility 
of a complete self-preserving solution with a K ~ t -1 asymptotic power law decay at 
high Reynolds numbers. In fact, this asymptotic decay law was found recently in the 
experiments of Walker (1986). 

(4) Since the assumption of complete self-preservation requires that G be constant 

— and since for high- Reynolds-number isotropic turbulence jg G > 2, whereas for low- 
Reynolds-number isotropic turbulence y$G <2 — it is clear that the entire process of 
isotropic decay from high Reynolds number initial conditions to the final period of decay 
cannot be described by the theory. The alternative conclusion drawn by George (1987, 
1989) arose from the use of a new scaling for the two-point triple velocity correlation that 
rendered a time dependent skewness where Sk ~ Rt ^ • However, this type of asymptotic 
behavior for the skewness wherein Sk — * 0 as — » oo is physically inconsistent since it 

yields zero energy transfer in the infinite Reynolds number limit. 

Within the framework of self-preservation, the physical origin of a K ~ t~ l power law 
decay becomes clear; it is the asymptotic state toward which a high-Reynolds number 
isotropic turbulence is driven in order to resolve an 0{^Ri ) imbalance between vortex 

stretching and viscous diffusion. The resolution of this imbalance also yields compatibility 
with Kolmogorov scaling. Results were presented which indicate that the complete self- 
preserving solution yields a better description of non-equilibrium isotropic turbulence 
than the commonly used turbulence models. It is also interesting to note that when 
the self-preserving assumption is extended to homogeneous shear flow, a production- 
equals-dissipation equilibrium can occur — preceded by a transient where K and e grow 
exponentially — as recently shown by Bernard and Speziale (1990). It thus appears that 
the theory of self-preservation in homogeneous turbulence has many interesting features 
that have not yet been fully understood and are worthy of further study. 
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Figure 1. Decay of turbulent kinetic energy in complete self-preserving isotropic turbu- 
lence for initial turbulence Reynolds numbers R to = 1000,5000 and 10,000: (a) C t2 — 
1.92, (b) C t2 = 5, and (c) C t2 = 8. 
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Figure 1. Decay of turbulent kinetic energy in complete self-preserving isotropic turbu- 
lence for initial turbulence Reynolds numbers R to = 1000,5000 and 10,000: (a) C i2 = 
1.92, (b) C (2 = 5, and (c) C t2 = 8. 
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Figure 1. Decay of turbulent kinetic energy in complete self- preserving isotropic turbu- 
lence for initial turbulence Reynolds numbers R to = 1000,5000 and 10,000: (a) C C2 = 
1.92, (b) C t2 = 5, and (c) C t2 = 8. 
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Figure 2. Decay of turbulent kinetic energy for a variety of initial conditions on G 

(C t 2 = y^G 0 ): self-preserving solution; asymptotic solution K ~ t~ a where 

a = (C e2 — l) -1 for C C2 < 2 and a = 1 for C C2 > 2. (a) = 1000, (b) R to = 10,000 

and (c) R to = 100,000. 
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Figure 2. Decay of turbulent kinetic energy for a variety of initial conditions on G 

(C e2 — y^Go): self-preserving solution; - - - asymptotic solution I< ~ t~ Q where 

a = (C (2 - I)" 1 for C C2 < 2 and a = 1 for C £2 > 2. (a) R to = 1000, (b) R to = 10,000 
and (c) Rt a = 100,000. 
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Figure 10. Decay of turbulent kinetic energy: o Self-preserving solution; 
(a) C t2 = 4, (b) C t2 = 4.5, and (c) C ej = 5. 
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